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ABSTRACT 

We discuss the applicability and derivation of window functions for cosmic microwave 
background experiments on large and intermediate angular scales. These window functions 
describe the response of the experiment to power in a particular mode of the fluctuation 
spectrum. We give general formulae, illustrated with specific examples, for the most com- 
mon observing strategies. 
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1. Introduction 



It has become conventional in cosmic microwave background (CMB) studies to de- 
scribe the sensitivity of experiments by "window" or "filter" functions. These functions 
describe the response of an experiment to the power in a particular mode of the underlying 
fluctuation spectrum. Plots of the window functions for various experiments are becom- 
ing common (see e.g. Bond 1990, Crittenden et al. 1993, Gorski 1993, White, Scott & 
Silk 1994), and are increasingly used to compare different experiments. We believe that 
it is useful to understand the derivation and meaning of these window functions from a 
generic point of view. Furthermore there are experiments for which analysis by a window 
function is too complicated to be useful, and it is important to understand what features 
of an experiment lead to this situation. 

We begin by assigning an "ideal" temperature T(n) to every point on the sky; this is 
the temperature that would be measured by a perfect experiment with an infinitely thin 
beam. Here n is a unit vector which can be expressed in the usual way in terms of the 
polar and azimuthal angles: 

n = (sin6>cos(/>, sin^sin^, cos6>). (1) 
We can expand the ideal temperature pattern in spherical harmonics: 

oo i 

T ( n ) = EE a>tmY tm (n) , (2) 

1=1 m=-£ 

where we have removed the constant (£ = 0) term. In any theory of fluctuations, the ae m 
are treated as random variables; if we assume nothing more than rotational invariance, we 
must have 

( a £m a £'m') ens = C l $m'm , (3) 

where the angle brackets denote an average over the statistical ensemble of temperature 
fluctuations. If the fluctuations are gaussian, all higher point autocorrelation functions are 
given in terms of the two-point function, and the set of C/s exhausts the content of the 
model. Thus, any gaussian model is completely specified by its predictions for the values 
of the CVs, which will depend on some underlying parameters. It is these parameters 
which we would like to measure with experimental data. The computation of the CVs has 
been discussed by several authors, including Peebles & Yu (1970); Wilson & Silk (1981); 
Bond & Efstathiou (1984,1987); Vittorio & Silk (1984,1992); Holtzman (1989); Sugiyama 
& Gouda (1992); Dodelson & Jubas (1993); and Stompor (1994). In terms of the C/s, the 
two-point autocorrelation function of the ideal temperatures T(n) is given by 

<T( ni )T(n 2 )) ens = — £(2£ + 1) Ci Ptim -n 2 ) , (4) 

i=\ 

where Pe(x) is a Legendre polynomial. 
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Of course no real experiment can measure the ideal temperatures. In fact, each ex- 
periment assigns a temperature, or a temperature "difference" , to points on the sky in a 
way which is generally unique to that experiment. Call the temperature or temperature 
difference assigned to a point n by a particular experiment T(n). In general it will be 
linearly related to the ideal temperatures in some neighborhood of n: 

T (n) = J dQ n/ M(n, n') T(n') , (5) 

where the mapping function M(n, n') depends on the detailed experimental strategy. The 
mapping function is usually too messy to compute in closed form, but it gives us a common 
way of thinking about experiments, and as we will see it is closely related to the more 
commonly used window function. 

There are a number of different aspects of the experiment that go into the mapping 
function. One common to all experiments is the beam profile function B(n,n'), which 
accounts for the directional response of the antenna. For large scale experiments, such 
as COBE (Smoot et al. 1992) and FIRS (Ganga et al. 1993), which map the sky this is 
the only effect, and the mapping function M(n, n') is equal to the beam profile function 
S(n,n'). 

However, experiments on smaller scales are usually more complicated. They typically 
use a chopping strategy, possibly coupled with a smooth scan. All cases are encompassed by 
the following treatment. We must specify the beam position function n(£), or equivalently 
9{t) and (f>(t), which tells us the position of the center of the beam at time t. We must also 
specify the weighting or lock-in function L(t), which tells us how different portions of the 
beam trajectory are weighted in computing an experimental temperature, and the overall 
normalization N. 

Now consider a particular time interval, labelled by i, which runs from t = U — to 
t = ti + 7}5i. To this time interval we assign an average position Tij, given in terms of 6i 
and §1 by eq. (1), where 

0i = Y dte(t), 

Oi Jti-Si/2 

(6) 

^ rti+Oi/2 
5i Jti-Si/2 

To the average position we assign the temperature 

jy pti+Si/2 
-Si/2 



1 rti+Si/2 
<t>i= T dt<f>(t). 
°i Jti-Si/2 

ssign the temperature 

~ N r t i+ s i/ 2 r 

T(n,) = - / dtL(t) / dtt n , B(n(t), n') T(n') . (7) 

Vi Jti-Si/2 J 



Eq. (7) is completely general, and applies to all experiments. It gives an implicit definition 
of the mapping function M(n, n'); compare eq. (7) with eq. (5). Below we will discuss a 
variety of possible choices for the beam position function n(t) and the lock-in function 
L(t). 
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To summarize, an experiment is completely specified by giving the beam profile func- 
tion -B(n, n'), the beam position function n(t), the lock-in function L(t), and the normal- 
ization N. All four must be given explicitly before the results of an experiment can be 
analyzed or understood. A measured value of "AT" without specification of all four of 
these experimental ingredients cannot be interpreted. 

Given the mapping function M(n, n'), we can compute the window function Wg(n, n') 
as follows. Let us begin with the two-point autocorrelation function of the experimental 
temperatures: 

(f ( ni )T(n 2 )) ens = J d n n[ J d^M^n^) M(n 2 ,n' 2 ) (T(n[)T(n' 2 )) ens 

= I dn n[ I rfO n ^M(n 1 ,n / 1 )M(n 2 ,n , 2 )^-f;(2£ + l)C,P,(n / 1 -n^) 

J J n 1=1 

1 oo 

= ^J](2^+l)Q^(ni,n 2 ), (8) 

where the last equation defines the window function: 

Wi(m,n 2 ) = J dn n[ J ^Mln^DM^n^^ni.^). (9) 

Often the window function is plotted in the literature as a function of I only. This case 
corresponds to the window function at zero-lag, W^(n, n). This is usually independent of 
the choice of n to a very good approximation. To specify slightly, we will assume that 
Wg(n, n) is indeed independent of n, and use the shorthand notation Wg = Wg(n, n). We 
will, however, also be interested in the complete window function Wg(iii, n 2 ). 

2. Simple Window Functions 

2.1 The beam profile 

As already noted, the simplest example of a window function is that which arises when 
the finite size of the beam is the only effect, as is the case for COBE (Smoot et al. 1992, 
Wright et al. 1994) and FIRS (Ganga et al. 1993). In this case, the mapping function 
M(n, n') is equal to the beam width function S(n, n'). If the beam profile is isotropic, 
then S(n, n') is a function of n-n' alone, and we specialize to this case from here on. We 
can expand S(n, n') in Legendre polynomials: 

B(n, n') = — J2( 2i + !) B l ^(n-n') , (10) 

71 £=0 

and rewrite the Pg, using the addition theorem for spherical harmonics, in eq. (9) to find 
that the window function is simply 

oo 

Wg(n,n') = J2B 2 gPi(n-n'). (11) 
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What we call B\ is called G e by Wright et al. (1994). For a gaussian beam profile, 



B(n,n') 



1 



exp[-9 2 /2a 2 ] , 



(12) 



27ra 2 



where 9 = cos 1 (n-n / ), we have to a very good approximation (Silk & Wilson 1980, Bond 
& Efstathiou 1984, White 1992) 



In general the effect of the finite beam width is to provide a high-£ cutoff at scales of the 
beam size I ~ a~ x . We note in passing that uncertainties in the value of a for a gaussian 
beam profile, or more generally the shape of the beam profile, can result in significant 
uncertainties in comparing experiment with theory, especially if the high-£ cutoff is in a 
range where Ci is changing rapidly with I. 

2.2 Constant elevation scans 

We now turn to small scale experiments which use nontrivial beam position functions 
n(t) and lock-in functions L(t). For these, it is possible to significantly simplify eqs. (7) 
and (9) only if the scan is performed at a constant 9. (Note that this needs to be the case 
only for one particular choice of coordinates, which need not be equivalent to any of the 
usual choices.) In this case the complete window function W^(n, n') is a function only of 
|0 — (f)'\, an enormous simplification. We specialize to this case for now, and will return to 
discuss the general case later. 

We first consider a stepped (as opposed to smooth) scan. In this case the beam is 
centered at a particular point n = (6>o, 0o) on the sky, and then "chopped" back and forth 
in the (^-direction. The instantaneous beam position n(t) is given by 



where «o is half of the peak-to-peak chop angle, and uj c /2tt is the chop frequency. (In 
practice it is a few Hertz, but the window function turns out to be independent of ui c .) 
Note that the angular separation on the sky is measured by 4>sm9o. 

In the case of a smooth scan, the beam is swept smoothly, at an angular velocity of 
cu s , in addition to being chopped. The instantaneous beam position n(t) is now given by 



and the data must be binned, as in eq. (7), by integrating t over the duration time 5 of 
a bin. It practice 5 is always a multiple of the period of the chop; that is, u c 5/2ti is an 
integer. 



B e (a) =exp[-±^+l) ( r 2 ] . 



(13) 



9{t) = 9 , 

4>{t) = 4> + a sin(w c t) , 



(14) 



9(t) = 9 , 

4>{t) = 4> + u s t + ao sin(o;ct) , 



(15) 
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We are now in a position to compute the window function, assuming either eq. (14) or 
eq. (15) for the instantaneous beam position. We further assume (again always the case in 
practice) that the lock-in function L(t) has the same periodicity as the chopping function, 
and that uj s <C uj c . Again making use of the addition theorem for spherical harmonics, we 
ultimately find from eqs. (7) and (9) that 

A 1 

WtW = N 2 B 2 (a) l*WM)| 2 L 2 m (a ) S 2 m (AcP) cos(m0) , (16) 

m=— I 

where, up to an irrelevant phase, 



L m (a ) = ^ [ +n/UJC dtL{t) e ~sin(u, c t) 

^ J-*/ Ut . 



and 

S „(A^i,W/2) = ^f (18) 

for a smooth scan. Here A0 = lv s 5 is the size of the bins in 0. For a stepped scan, S m = 1. 
For more details in the context of specific choices of L(t), see Dodelson & Jubas (1993); 
White et al. (1993). 

If we can neglect the curvature of the line and assume that it is an arc of a great circle 
(usually a very good approximation), then we can set 9 = 7r/2 and let 4> be the angle on 
the sky; in this case 



(2£-2r)!(2r) 



WW) = N 2 B 2 (a) £J V ^ r , ( 7^,'p Lt 2r (a ) Sl 2r (A<P) cos[(£ - 2r)0] (19) 

which is easy to implement numerically. Note that now 0, A0, and cto are all defined as 
angles on the sky. 

2.3 The lock-in 

The simplest lock-in function is that for a "square wave chop" recently used by Tener- 
ife (Hancock et al. 1994), MSAM (Cheng et al. 1994), OVRO (Myers, Readhead & 
Lawrence 1993) and Python (Dragovan et al. 1994). In this strategy the temperature 
assigned to 4>q is a weighted sum of temperatures along a line, which we assume to be of 
constant elevation. The telescope moves rapidly between the observed points, stopping 
and taking data at set positions on the line. The weights assigned to points on the sky for 
three different "switching strategies" are 



(0 - 4>o)/do 


-1 


1 
3 





+1 


+i 


2-beam 


+1 








-i 


3-beam 


l 

2 




+1 




i 

2 


4-beam 


+1 


3 
4 




+! 


1 

4 
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In our notation such a strategy is implemented by taking L(t) to be a linear combination 
of Dirac delta functions. For a two-beam, three-beam, or four-beam switching strategy, 
we have 



{+ 6(t - t c ) - S(t + t c ) (2-beam) 
-h6(t-t c ) + 6(t)-±6(t + t c ) (3-beam) (20) 

+ W - tc) - \&it - # c ) + lS(t + # c ) - W + t c ) (4-beam) 

where £ = sin -1 (1/3) and t c = n/2u c is the time to chop from = 0o to 4> = 4>o + "o- 
In this case the mapping function M(n, n') defined through eq. (7) reduces to a weighted 
sum of beam profile functions -B(n, n'). From eq. (17) we find immediately that, up to an 
irrelevant phase, 

2 sin(m«o) (2-beam) 
L m (a ) = { 2sin 2 (ma /2) (3-beam) (21) 

|[sin(mao) — 3 sin(|mo;o)] (4-beam) 

Notice that L m (a ) scales as Oq -1 for the n-beam switching strategy. In general the 
window function for any kind of differencing experiment is suppressed at low £, since any 
long wavelength perturbation is removed by the differencing. Since the low t cutoff is 
controlled by ao while the high I cutoff is specified by a, one can increase both the height 
and width of Wg by separating these scales as much as possible. 

To make contact with forms of Wg frequently quoted in the literature, we note that we 
can substitute eq. (21) into eq. (19) with = (zero lag) and use the addition theorem for 
spherical harmonics to obtain 

2 [1 - P^(cos 2oq)] (2-beam) 
W £ = Bj{a) { \ [3 - 4P £ (cosa ) + P^(cos2a )] (3-beam) (22) 

| [10 — 15P£(cos |a ) + 6P^(cos |a ) ~ -P^cos 2a )] (4-beam) 

Three other illustrative choices of the lock-in function for differencing experiments are 
the "square-wave lock-in," 

Lit) = 2 sgn[a sm(w c t)] , 
L m (a ) = 2H Q (ma Q ) , 

used by SP91 (Gaier et al. 1991) & ARGO (de Bernardis et al. 1994), the "sine- wave 
lock-in," 

L(t) = 7ism(u c t) , ^ 
L m (a ) = ir Ji(ma ) , 



used by MAX (Gundersen et al. 1993, Meinhold et al. 1993) and the double-angle "cosine 
lock-in" , 

Lit) = tv cos(2uj c t) , 
W V / (25) 

i m («o) = 7r J 2 (ma ) , 

used by Saskatoon (Wollack et al. 1994). Here Hq(x) is the Struve function, and J n {x) are 
Bessel functions of the first kind. The numerical prefactors are chosen so that 



2ttJ_ 



+ 7r/uj c 

dt\L(t)\ = 2, (26) 

t/lj c 



which is a common way of normalizing an experiment (more on this below). Sometimes 
the L m (o!o)'s of eqs. (23) and (24) are approximated by the L m (ao) for a "2-beam chop", 
as given in eq. (20), but this can result in significant errors. For example, after taking 
into account the normalizations used by these experiments, we find that for SP91 the 
approximation is off by 20% (see Dodelson & Stebbins 1994) while for MAX it differs by 
10%. Also, the Saskatoon lock-in function of eq. (25) is only roughly approximated by the 
3-beam result of eq. (21). 

2.4 The normalization 

Finally we consider the normalization factor N. For SP91 and MAX this is chosen so 
that if there is a sharp boundary between two regions of constant temperatures 7\ and 
T2, then aiming at a point on the boundary, and chopping perpendicular to the boundary, 
gives T = T 2 — T\. For SP91, the normalization is computed assuming a perfect, point-like 
beam, corresponding to a = in eq. (13). For a perfect beam, the lock-in factor L m (ao) 
given in eqs. (21-25) is already normalized so that N = 1. However, if the normalization 
is done assuming the actual beam profile of the experiment (as is the case for MAX and 
Saskatoon), then there are corrections which must be computed. The result for MAX was 
presented in Srednicki et al. (1993), but here we give a more general treatment. 

We first make the "flat sky approximation" near the T\-T 2 boundary, which we shall 
take to be the line of longitude <fi = 0. We treat x = <j> and y = ^ — 9 as cartesian 
coordinates. For n = (0, 0) and n' = (x, y), eq. (12) for the beam profile becomes 

B(x,y) = ^e- (x2+y2)/2a2 ■ (27) 

Eq. (5) can now be written as 

/+00 
dxdyM(x,y)T(x,y), (28) 
-00 

where the mapping function is given by 

Mi \ NuJc [ +n/ " C * r fA ( [x-aosm(u; c t)} 2 + y 2 \ 
M (x,y) = /n , 9 o / dtL(t)exp( — . (29) 
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Take the temperature profile to be T(x, y) = T 9(x), where 9(x) is the step function, and 
demand that T(0, 0) = T . For the MAX lock-in function given by eq. (24), we get 



/ + 7T 
dr sin r erf (7 sin r) 
(30) 

= N(^^)iFi(i 2; - 7 2 ), 

where 7 = cto/v^ ", er f is the error function, and 1F1 is the confluent hypergeometric 
function. For MAX, with a = 0°65 and a = 0.425 x 0°5, this gives A 2 = 1.13. (Note 
that it is A 2 which appears in the window function.) 

For Saskatoon, the only change is that the normalizing temperature profile is Tq where 
M(x, y) > 0, and zero elsewhere. The mapping function follows from substituting eq. (25) 
into eq. (29), which also implicitly defines the region where T(x, y) 7^ in eq. (28). Setting 
f (0, 0) to T as before gives 

N f , f* , / n \ ( [x-a smr} 2 \ 
1 = = / dx / drcos(2r)exp -± — - . (31) 



2aV2^ J J-n V 2d 



where the integration is over the region of x where M(x, y) > 0. With a = 2°45 and 
a = 0.425 x 1°44, we find A 2 = 1.74, quite a large correction. Also note that the mapping 
M(x, y) defined by eq. (25) and eq. (29) is not well approximated by a square wave chop 
pattern, which would be three gaussian beam profiles at x = — ocq, and cto with weights 
— i, +1 and — I respectively. 



3. Other Strategies 



One other "differencing" strategy that has been proposed recently is used by the "White 
Dish" experiment (Tucker et al. 1993), which assigns to the point n a temperature T(n) 
which is given by a particular weighted average of measured temperatures in a circle around 
that point. The actual strategy used is very difficult to model effectively. However, if we 
modify their "Method II" analysis to neglect binning, then it is straightforward to compute 
Wi = Wg(n, n). While the off-diagonal elements cannot be simply constructed, a numerical 
procedure similar to that used in Srednicki et al. (1993) would be feasible. 

To get Wi, we first rotate coordinates so that n is at 9 = with the circle being in <fi at 
fixed 9 — 9q. The analogue of eq. (7) is now to extract the nth harmonic of the temperature 
around the circle, which corresponds to a window function of the form [c.f. eq. (16)] 



(32) 
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where in the last line the limit ^ <C 1 has been used. The term in parenthesis on the 
last line is very close to unity for I ^> n. Assuming that a temperature profile which 
is Tq for both < (p < n/2 and 7r < <p < 37r/2 and zero elsewhere would be assigned a 
temperature difference of T = To, when n = 2, we get N = 1. A less artificial normalization 
would require more information than is specified in the paper (Tucker et al. 1993). The 
White Dish experiment has 9q = 14', a = 0.425 x 12', and n = 2, with the resulting 
"temperatures" binned into 4 positions T lr ..,T 4 in a square of side 6' = 23.6'. The four 
temperatures are assigned to consecutive corners clockwise around the square (Tucker et 
al 1993). 

Eq. (32) can be compared with the window function for a "square wave chop" procedure 
where we simply sum the temperatures in a square: T = \{Ti — T 2 + T 3 — T 4 ), neglecting 
how they were measured, to obtain 

1 - 2P £ (cos 6') + P e (cos y/29') 

Assuming that N = 1, the two window functions peak in the same place (£ ~ 500), but 
differ by 20% at the peak and have a different scaling with £ off the peak. Additionally, 
neither of these methods accurately reflects the correlations induced by coarse binning of 
the data. 

With an analysis procedure as difficult to model as White Dish, the window function 
approach is of limited utility and one should resort to Monte-Carlo simulations of the 
observing strategy, for each theory being tested. Alternatively the applicability of a window 
functions should be kept in mind when the analysis procedure is designed. 

4. Scans at Varying Elevation 

If data points are not taken at constant elevation, such as in the GUM scan of MAX 
(Gundersen et al. 1993), and chopping is used, computing the window function at nonzero 
lag (n ^ n') is impossible analytically. The information which is needed to compare data 
with a theory is the autocorrelation function of the experimental temperatures, eq. (8), 
computed with the C/s of the theory in question. It is then usually easier to compute 
this directly, using numerical methods, than to compute the complete window function 
Wi(n, n'). However, even numerical methods become cumbersome if the observing strategy 
is complex. Accounting for chopping generally requires that a double integral be done 
numerically (Srednicki et al. 1993). If there is, in addition, a smooth (as opposed to 
stepped) scan, then a quadruple integral must be done numerically, and this is not feasible 
in general. If the data is binned finely enough, then the effects of the smooth scan are 
small, and this is not a problem. Here "finely enough" means A0 < a , fortunately the 
case for MAX-GUM. 

5. Conclusions 

Using the formulae presented in this letter it is possible to compute window func- 
tions for most of the current large and intermediate scale CMB anisotropy experiments. 



W £ = B 2 e (a) 
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We stress however that there exist several generic cases in which the window function 
approach is not the optimal method of analysis. These are when the correlation matrix 
(T(ni)T(ii2)) is anisotropic (such as in the case of the GUM scan of the MAX exper- 
iment or multiple scans of the SP91 experiment) or if the experimental procedure makes 
analytic calculation of the window function difficult. In these cases, it is generally much 
easier to compute eq. (8) directly (by numerical methods), or to do a Monte Carlo analysis 
of the experiment, than it is to try to calculate the window function We(n, n'). 

Both diagonal and off-diagonal window functions are easy to construct for experiments 
in which the scanning and chopping directions are constant, regardless of the type of chop- 
ping (square wave, sine, cosine), and regardless of whether the scan is smooth or stepped. 
For these cases, the window function provides an efficient method for comparing theory 
with data. However, with the current refined state of CMB anisotropy measurements, it 
is important to used the right window function. A simple stepped-scan, square-wave chop 
approximation to all experiments is no longer accurate enough for the quality of data now 
available. 

We would like to thank Douglas Scott for useful conversations and comments on the 
manuscript. This work was supported in part by NSF Grant Nos. PHY-91-16964 and 
AST-91-20005. MW. acknowledges the support of a fellowship from the TNRLC. 
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